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Abstract. Published survival data for V79 cells irradiated by monoenergetic protons, 
helium-3, carbon, and oxygen ions and for CHO cells irradiated by carbon ions 
have been analyzed using the probabilistic two-stage model of cell inactivation. 
Three different classes of DNA damages formed by traversing particles have been 
distinguished, namely severe single-track damages which might lead to cell inactivation 
directly, less severe damages where cell inactivation is caused by their combinations, 
and damages of negligible severity that can be repaired easily. Probabilities of single 
ions to form these damages have been assessed in dependence on their linear energy 
transfer (LET) values. 

Damage induction probabilities increase with atomic number and LET. While 
combined damages play crucial role at lower LET values, single-track damages 
dominate in high-LET regions. The yields of single-track lethal damages for protons 
have been compared with the Monte Carlo estimates of complex DNA lesions, 
indicating that lethal events correlate well with complex DNA double-strand breaks. 
The decrease in the single-track damage probability for protons of LET above approx. 
30 keV//xm, suggested by limited experimental evidence, is discussed, together with 
the consequent differences in the mechanisms of biological effects between protons and 
heavier ions. Applications of the results in hadrontherapy treatment planning are 
outlined. 
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Hadron radiotherapy, based on irradiating tumours by beams of accelerated protons 
and other ions, is e xpected to significantly increase the cure rate of cancer in near 
future; compare e.g . ( Maver et al 112004 , Krengli and Orecchiall2004 , Brahme et al IhoOl 



Debus et al 1998). The well-known rationale for using proton beams in radiotherapy 



lies in the characteristic pattern of their energy deposition to matter, the Bra gg peak. 



leadi n g to the possibili t y of achieving hig hly conformal dose distributions (W ilson 



1946, 



Zurlo et al 



2000, 



Cella et al 1 120011 ). Beams of light ions possess additional 



advantages: They are characterized by higher linear energy transfer (LET) values, 
leading to enhanced relative bi ological effectiveness (R BE) and diminishing oxygen 
enhancement ratio (OER), e.g. ((Wambersie et al 1 120041 ) . They are therefore expected 
to be especial ly suitable for curing t umours resistant to conve ntional radiotherapy; 
compare e.g. ([Wambersie et al 2004 . Schulz-Ertner et al 2002 ). Another advantage 
consist s in the possibility of online dose monitoring by positron emission tomography, 
PET ((Enghardt et al I 120041 ) . Furthermore, recent clinical studies indicate that 



unconventional fractionation schemes might b e used for selected tumours, reducing the 
overall treatment time to a few sessions only (jTsujii et al 1120041 ) . 

Several ded icated hadrontherapy facilities have been launched recently, compare 
(Sis tersonl 12005) . In most of them only proton beams are available; however, in several 
centres, existing or being built, protons as well as light ions up to carbon or oxygen may 
be used. For the choice of an optimal ion in given clinical situations, analyses of energy 
loss, scattering and fragmentatio n phenomena for different ions and comparisons of their 



dose distributions are necessary ((Brahme et al 



2001 



Kramer et al 



2000, 



Schall et al 



1996). These physical analyses have to be complemented then by biology-oriented 



studies, analyzing the mechanisms of biological effects for different ions using available 
data and adequate models. 

However, existing hadrontherapy treatment planning approaches have been based 
on detailed description of the physical processes only, and have not addressed the 
biological phase of the underlying radiobiological mechanism in detail. E.g. the 
treatment planning procedure used at HIMAC (Chiba, Japan) is based on interpolating 
the a, (3 coefficients of the linear-quadratic (LQ) fits to experimental survival data, and 
on the similarity in biolog i cal eff ects between car bon and neutron bea ms of similar LET 
values (jKanai et al II1997L Il999l ): compare also (jKagawa et al 1 120021 ) . In the approach 
used in clinical hadrontherapy applications at GSI (Darmstadt, Germany), the local 
effect model (LEM), the effect of ion tracks on a microscopic scale is assumed to be equ al 
to that of correspondingly high photon doses (compare ((Kramer and Schola 120001 ) ) : 
again, the underlying biological processes of damage induction and repair have not been 
taken into account explicitly. On the contrary, the present paper summarizes the results 
of a detailed biology-oriented analysis of cell inactivation data for proton s, helium-3 
carbon and oxygen ions . Experimental data gathered by several g r oups (Belli e t al 



1998. 



Perris et al 



1986, 



Goodhead et al 



199; 



Folkard et al 



1996, 



Wevrather et al 
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1999L IStoll et al Ill995) have been analyse d with the help of a detailed probabilistic 



radiobiological model (|Kundrat et al 1120051 ). 



2. Probabilistic two-stage model of radiobiological effects 

The processes running in a cell after irradiation (under conditions usual in radiotherapy) 
proceed in two distinct, subsequent phases. The first stage includes processes running 
immediately after the impact of individual ionizing particles, i.e. energy transfer events, 
radical formation and diffusion, chemical reactions and formation of DNA damages. 
The other stage is formed then by the subsequent cellular response which follows as a 
reaction to the effect of all ionizing particles having hit a cell nucleus or chromosomal 
system at a given dose; this stage includes the processes of damage repair or misrepair 
and further biological processes leading to cell survival or inactivation. Distinguishing 
these two phases represents the basis of the probabilistic two-stage model of biological 
effects of ionizing particles. 

In the full scheme of the model, the following characteristics are taken into account: 
(i) the stochastic distribution of particle tracks over the irradiated cell population, (ii) 
the distribution of transferred energy, (hi) the probability of single traversing particles 
to induce DNA damages of different severity, (iv) saturation or synergetic combinations 
of individual d amages, and (v) the e ffects of cellular repair systems. The model scheme, 



presented in (IKundrat et al 1 120051 ) . allows to describe different classes of survival 
curves, including also low-dose hypersensitivity phenomena, on the basis of detailed 
characteristics of the mentioned processes, especially the effects of damage induction 
and repair processes. However, in the following comparison of the effects of different 
light ions that are expected to be used in hadrontherapy, a simplified model scheme will 
be used, taking into account only DNA damages that are practically unrepairable. 
The cell inactivation probability at applied dose D is given by 

s(D) = Y / Pk(D)q k . (1) 

k 

Here, the number k (k > 0) of primary particles (particle tracks) traversing individual 
cell nuclei is given by Poisson distribution 

P k {D)= i ^eM-hD) , (2) 

the average number per 1 Gy, h, being proportional to the geometric cross-section of 
the nuclei, a. The survival probability after k particles have traversed the nucleus has 
been denoted by It can be calculated by 

Qk = 1 - E({)a'(l - a)" " b) k ~^[l - r#] , (3) 

i=0 j=0 

where combinatorial numbers (*) = juj^Zty ■ Here a and b stand for the induction 
probability of different damages and rf- for repair probabilities of their different 
combinations; Tqq = 1 = (compare below). Three different classes of damage formed 
by individual particles have been distinguished: 
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(i) Severe damage formed by a single track that is capable of causing cell inactivation 
even if this is the only damage to the cell ("single-track" or "single-particle induced" 
damage). The probability of inducing such a damage, per track, has been denoted 
by a. 

(ii) Damages of lower severity, at least two of which must combine to be lethal 
("combined" or "two-track" damages). The probability of their formation (per 
track) has been denoted by b. 

(iii) Lesions that do not represent any significant challenge to the cell and its repair 
systems. Such negligible lesions have not been included in Eq. ©. 

The above classification of DNA damage is illustrated in Figure H 



negligible 



lower severity (b) 



severe (a) 



► 100% 

Figure 1. Schematic illustration of the damage classification according to damage 
complexity. 

Eq. Q represents the general formula for cell inactivation after a passage of k 
particles. The effects of cell repair systems have been included in repair probabilities rf-. 
In the present analysis, however, only unrepairable damages have been considered, for 
which rfj = 0. In other words, the probabilities a and b derived here represent the lowest 
estimates of damage induction at corresponding energy transfers (LET values), since 
unrepairable damages form, essentially, subsets of all (both repairable and unrepairable) 
damages in the above classification (Eq. (jHJ), Figure HJ). 

Considering unrepairable damages only, cell inactivation follows, then, if at least 
one single-track or at least two combined-type damages have been formed by the k 
particles traversing the nucleus. Indeed, for rf- = 0, Eq. (JHJ) may be simplified to 

q k = (l- a ) k -(l-a) k ^2())b>(l-b) k - j = (l-a) fc [(l-6) fe + A;6(l-6) fc - 1 ](4) 

(with go — 1 and qi — 1 — a), i.e., the cell survives if no single-track and not more than 
a single 6-type damage have been formed by the k particles. $ 

| The corresponding Eqs. (9-f5) in I Kundrat et al 112005^ represented simplified versions of formulas 



derived from Eq. as combined damages b were taken into account in a simplified manner 

only: The synergetic effects were represented by terms (1 — fr 2 )*-'^- 1 )/ 2 only, yielding e.g. = 
(1 — a) k (l — 6 2 ) fe ( fc ~ 1 )/ 2 after neglecting the repair probabilities. By comparing this formula with 
Eq. 0} it is clear that the role of combined damages was slightly underestimated, especially in the 
region of higher particle numbers. I.e. the analyses based on precise formulas derived from Eq. © 
yield slightly hig her combined -damage induction probabilities as compared to the simplified formulas 
used in |Kundrat e.t al I Eof).^ . This fact, however, does not affect the conclusions drawn in that 
paper, especially those concerning the role of single-particle and combined damage induction and repair 
processes with respect to the shapes of survival curves. 
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Table 1. Physical parameters of particles in the analyzed data sets. Energy, 
LET values and ra nges as given b y the authors , or calculated using PSTAR tables 
i Bereer et afl EoOO) or SPJM-2003 jZierierlkooij) . 





Energy 


LET 


Range in water 


Reference 




[MeV/u] 


[keV/^m] 


[mm] 




p 


7 1 
1 .4 


O.O 


n 70 


rcrris et al iiyoo) 




O.U1 


7 7 


n Q7 
v. 61 


rjcin et al I laao I 




3.66 


10.1 


0.21 


Folkard et al ( 1996) 




3.20 


11.0 


0.17 


Belli et al (1998) 




3.0 


11.7 


0.15 


Perris et al ( 1986) 




1.83 


17.8 


0.065 


Folkard et al ( 1996) 




1.41 


20.0 


0.043 


Belli et al C1998) 




1.4 


on q 


U.U4Z 


v^ooancacL et al iiyyz) 




1.16 


23.0 


0.040 


Goodhead et al (1992) 




1.07 


27.6 


0.027 


Folkarii <pf a( (1996) 




0.76 


30.5 


0.016 


Belli et (1998) 




f\ P A 

0.64 


34.6 


a m o 

0.012 


Bclh^^g^ ( 1998 ) 




0.57 


37.8 


0.010 


Belli (1998) 


3 He 


2.30 


58.9 


0.074 


Folkard et al (1996) 




1.39 


88.3 


0.033 






1.13 


105.8 


A AO A 

0.024 




C 


266.4 


13.7 


140 


Wevrather et al (1999) 




190.7 


16.8 


80 






76.9 


32.4 


15. 9 






18.0 


103 


1.2 






11.0 


153.5 


0.5 






5.4 


275.1 


0.15 






4.2 


339.1 


0.1 






2.4 


482.7 


0.045 







396 


18 


202 


Stoll et al (1995) 




88 


46 


15.2 






10.7 


238 


0.36 






1.9 


754 


0.030 





3. Analysis of experimental data 



Survival data for V79 cells irradiated by protons at 0.57 
Folkard et ffl/~lll996 



flFolkard et al 



Penis alUlQS O 



7.4 MeV (jBelli al 111998 . 
Goodhead fiiT^71ll992h. helium-3 at 3 .4 - 6.9 MeV 



1996|), ca rbon at 2.4 - 266 MeV/u (|Wevrather et a/ 1119 99) and oxygen 
ions at 1.9 - 396 MeV/u (jStoll et a l 1995) have been analysed using the model reviewed 
in the preceding section. The list of particle energies, LET values and ranges in water 
is given in Table H 

To facilitate a systematic analysis of the given data sets, the damage induction 
probabilities for a given ion have been considered in dependence on LET value A 
instead of deriving their values for individual measured survival curves independently. 
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These LET dependences have been represented by flexible test functions involving a low 
number of auxiliary parameters, a _2 and &o-2 : 

a(X) = a„(l - exp(-(a 1 A) a2 )) , b(\) = b {l - exp(-(6 1 A) fe )) . (5) 

Parameter values have been determined using standard optimization methods; mainly 
the SIMPLEX an d MIGRAD methods implemented in the MINUIT multivariate 
minimization tool (j.Tameslll994[ ) have been applied to dedicated computer codes written 
in FORTRAN. Weighted least-square method has been use d to construct the objective 
function, x 2 , as described previously ( Kundrat et al 200.4 ). 

Model representations of survival curves for V79 cells irradiated by different ions 
are shown in Figures EE In Figure H3 the probabilities of single-track and combined 
damage induction are plotted as functions of LET for different ions. Values of auxiliary 
model parameters are listed in Table El 





Folkard 10 keV/micrometer - 

Belli 1 1 keV/micrometer 
Perris 1 1 .7 keV/micrometer - 



D[Gy] 





Belli 34.6 keV/micrometer - 
Belli 37.8 keV/micrometer 



1 2 3 4 5 



Figu re 2. Survival curves for V79 cells irradiated by protons; data taken from 
( Belli et al Ill998l iFolkard et al lll996l IPerris e^"oF1ll986l iGoodhead et al Ill992 ). 






Figure 3. V79 cells irradiated by low-energy 3 He ions; data from IjFolkard et 
1996). 



A similar analysis has been performed for s urvival data of CHO-K 1 cells 
irradiated by carbon ions, measured experimentally bv IWevrather et al I ( 19991 ) . Model 
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13.7 keV/micrometer ■ 

16.8 keV/micrometer □ 



3 4 5 
D[Gy] 




1 2 3 



D[Gy] 



275.1 keV/micrometer ■ 
339.1 keV/micrometer □ 
482.7 keV/micrometer * 




2 3 4 5 6 7 

D[Gy] 



Figure 4. Representation of surv ival curves for V79 cells irradiated by carbon ions; 
data from I Wevrather et al 1 1999|) . 



18 keV/micrometer 
46 keV/micrometer 



238 keV/micrometer 
754 keV/micrometer 



1 2 3 4 5 
D[Gy] 




1 2 3 4 5 
D[Gy] 



Figure 5. Survival of V79 cells irradiated by oxygen ions; data measured bv lStoll et 
(1993). 
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Figure 6. Probabilities of single-track and combined-damage induction for protons 
(top left), helium-3 (top right), carbon (bottom left) and oxygen ions (bottom right) in 
dependence on LET. For helium-3, the combined damages might have been neglected 
at all, 6 = 0. 



representation of survival curves and the derived damage induction probabilities are 
shown in Figure the parameter values are given in Table |21 

The present results correspond to the geometrical cross section for V79 cells 
chosen uniquely for all data sets to be a = 87.8 jim 2 , based on the value reported 



Cell inactivation by protons and light ions 



8 



Table 2. Values of model parameters, describing the dependence of damage induction 
probabilities on ion LET values for V79 and CHO-K1 cells, Eq. JSJ. 





V79 


CHO-K1 




P 


3 He 


C 


O 


C 


a 




0.045 


0.25 


0.41 


0.80 


0.50 




[/im/keV] 


0.036 


0.015 


0.0063 


0.0042 


0.0064 


«2 




3.57 


2.79 


3.21 


4.17 


2.0 


b 




0.017 


0.0 


1.0 


0.39 


0.26 


h 


[/xm/keV] 


0.13 




0.0012 


0.0093 


0.0011 


b 2 




2.81 




1.0 


1.44 


0.69 




3 4 5 
D[Gy] 



3 4 
D[Gy] 



E 0.2 




50 100 150 200 250 300 350 400 450 500 
LET [keV/micrometer] 



Figure 7. Survival curves (data from i Wevrather et ^119 99^ and damage induction 
probabilities for CHO-K1 cells irradiated by carbon ions. 



in (jWevrather et al I I19991 ). in order to enable comparisons between different ions. 
Values in the range of 50 - 135 urn 2 , a lso reported in the literature (compare e.g. 
(Sc holz and Kraftl Il995l iBelli et al I Il998l )). have yielded similar results (figures not 
shown here). For CHO cells, the value of geometrical cross sectio n a = 108 /xm 2 was 
taken, as reported in the given experiment (jWevrather et al II19991 ). 
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4. Discussion 

4-1- Single-track damage induced by protons: Possible indications for differences in 
radiobiological mechanisms between protons and ions 



Meas ured data for protons of LET values above approximately 30 keV / fim ((Belli et al 



19981 ) suggest that the single-track damage formation probability, a, might decrease in 
this region; compare Figure |H1 where more detailed fits to the data are shown, based on 
the parameterization 

o(A) = a (l - exp(-(a 1 A) fl2 )) exp(-(a 3 A) a4 ) , 

b(X) =& (l-exp(-(M) 62 )) • (6) 



Where as the value of the goodness-of-fits criterion is x — 350 for the data of Bell i et al 



(1998) in the joint fits to proton- induced survival shown in Figure El the corresponding 
value for the detailed fits shown in Figure |H1 is Xdetaiied = 55. The goodness-of-fits to 
this data set can be, in fact, further improved if, in addition to a decreasing single- 
track damage probability, the combined-damages are allowed to domina te in the whole 



proto ns' LET range and their repair is taken into account, compare ((Kundrat et al 



2005J). 

To verify the hypothesis of decreasing single-track damage probability for protons of 
high LET, further experimental studies are necessary, as the existing evidence is limited 
to two survival curves measured under a single experimental setup only. Even though 
this decrease corresponds to very limited proton ranges in tissues only (approximately 15 
fim or smaller) and will probably not manifest in clinical applications, this finding would 
be important for understanding the mechanisms underlying the radiobiological effects: 
For heavier ions, no such complex behaviour has been observed in the studied LET 
ranges, the damage induction probabilities exhibiting saturation characteristics only 
(Figure EI). It means that for heavier ions the RBE effects termed often as "overkill" 
occur as a result of too high energy deposits to individual cells. For protons, on the 
other hand, this evidence suggests that they are already the effects of single particles 
on the level of DNA damage formation that get saturated and are responsib le for the 



decrease of biological effectiveness observed experimentally ((Belli et al 1119981 ) . 



4-2. The influence of repair processes 

As mentioned in the preceding paragraph, the effects of cellular repair processes might 
be important in representing the outcome of proton irradiation. Also the fits to 3 He 
data in the high-dose region might be improved significantly by incorporating the 
repair probabilities (Figure OJ). In this region, however, the experimental conditions 
are very challenging and, consequently, the data in this region have a limited reliability 



(Fol kard et a/lll996l ): therefore, the authors have interpreted the upward-bending of 
survival curves as artificial effects only. However, the results shown in Figure El 
demonstrate that similar upward bending might follow as a result of repair processes; 




10 15 20 25 30 35 40 
LET [keV/micrometer] 



Figure 8. Detailed fits to proton-induced V79 inactivation, indicating a decrease in 
the single-track damage for mation for protons at LET values above approximately 30 



keV//xm; data measured bv iBelli et al I |l998), calculations based on Eq. JHJl. 



compare also the general discussion in (jKundrat et al II2005T ). The damage induction 
probabilities are almost identical with those presented in Figure El (in both cases, the 
combined-damage probabilities have been neglected, 6 = 0); the survival values at high 
doses being mainly influenced. Again, additional experimental data would be helpful in 
solving this issue. 




3 keV/micrometer - 



D[Gy] 



D[Gy] 





20 40 60 80 100 
LET [keV/micrometer] 



LET 80 
[keV/micrometer] 



15 

k 



Figure 9. Upper panel: Representation of su rvival curves for V79 cells irradiated by 
helium- 3 ions (data from jFolkard et al Ill996|) 1 when repair r a of single-track damage 
taken into account. Lower panel: Damage induction and repair probabilities. 



The goodness-of-fits for heavier ions does not increase significantly when the repair 
processes are taken into account. For the sake of simplicity, only unrepairable damages 
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have, therefore, been considered in the present work (Section |3J). The induction 
probabilities for both single-track and combined damages have been estimated on 
the basis of survival data only. This slightly limits the reliability of estimating the 
roles of combined-track events, b, in the higher-LET regions, whereas the unrepairable 
single-track damages, a, have been assessed reliably in the whole studied LET regions. 
Somewhat higher values of damage induction probabilities would be obtained if the 
influence of repair processes were included. A corresponding concept, including a wider 
class of damages and their repair by the cells, has been used in analyzing radiobiological 
data for cell lines of different radiosensitivity ; the results will be presented and discussed 



in detail elsewhere (jHromcikova et al 1 120051 ) 



4-3. On the interpretation of the damage classification 

Clustered DNA lesions of high complexity (locally multiply damaged sites, complex 
lesions) are thought to represent significant challenges to cellular repai r systems, and to 



play crucial roles in cel l inactivation by ionizing particles (IWard lll985LlGoodheadl ll994. 



Ottolenghi et al 1997). The existing estimates of their yields have been based mostly 
on Monte Carlo calculations involving track structure simulations and simplified models 
of DNA and/or chromatin structure. Due to their complexity and computational costs, 
these approaches have been performed for photons and low-energy particles only. On the 
other hand, the present approach, based on analyzing survival data, enables to derive 
the yields of lethal events also for high-Z high-energy ions (FigurelBl). By comparing the 
present results with the mentioned Monte Carlo calculations, mechanistic interpretation 
of the damage classification used in the present work might be sought, and the fraction 
of lethal events among different classes of complex lesions can be estimated. 

As the first step in this direction, the present results for lethal damage induced 
by protons in V79 cells have been compared to estimates of complex lesions 
calculated on the basis of a fa st Monte Carlo damage simulation model proposed by 



Semenenko and Stewart I ([20041 ). In Figure EH the LET-dependent yields of single-track 
lethal events, a, are compared to the total number of double-strand breaks (DSBs) and 
DSBs composed of at least 4, 6, 8, and 10 elementary lesions (single-strand breaks 
(SSBs), base damages or abasic sites). This figure indicates that the present results 
are consistent with the trends predicted by Monte Carlo calculations. In Figure the 
ratio between the number of lethal events and DSBs of different complexity is shown, 
again as a function of LET. Whereas only 0.05 - 0.5% of all DSBs formed (denoted 
by DSB total in Figures ITUl and fTT[) are lethal, this portion increases to approximately 
3 - 4% for DSBs composed of at least 8 elementary lesions (DSB 8+); this number being 
uniform (LET-independent), i.e. being consistent with the hypothesis of a uniform repair 
probability for these specific damages. 

These results indicate that the single-track lethal damages, a, might be related to 
clustered lesions of very high complexity, and correlate best with DSB 8+, containing 
at least 6 elementary DNA lesions in addition to 2 SSBs forming the DSB. Similarly, 
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Figure 10. Comparisons of damage yields for protons of different LET values: single- 
track lethal damage as derived in the present work, and Monte Carlo estimates of total 
DSB numbers and yields of DSBs of higher complexity, consisting of at least 4, 6, 8, 
or 10 elementary DNA lesions (DBS 4+, 6+, 8+, 10+, respectively). 
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Figure 11. The ratio of lethal events among all DSBs (DSB total) and complex DSBs 
formed by at least 4, 6, 8, or 10 elementary lesions. 



the combined damages, b, might correspond to pairwise combinations of DSBs formed 
by different tracks, lead ing to large-scale chromosomal aberrations; compare e.g. 
(jSachs and Brenner! I1993T ). Further damage induction and repair studies are, however, 
necessary to identify the biophysical nature of lethal events reliably. Furthermore, this 
comparison has been performed for protons only, given the limitations of existing Monte 
Carlo approaches. Efforts will therefore be made to extend these calculations to heavier 
ions and higher energies, using correspondingly simplified Monte Carlo-based models. 



4-4- Applications in hadrontherapy treatment planning 



The analyzed proton and helium-3 d ata ([Belli et al I Il998l iFolkard et al 



1996 



Perris et al I ll98fiL lOoodhead et al I Il992h concern irradiation by low-energy particles 
only, with ranges in tissue of approx. 10 - 0.7 mm for protons and 24 - 74 /im for 
3 He ions, respectively (Table HJ), corresponding to track ends of these particles. On the 



Cell inactivation by protons and light ions 



13 



other hand, data for carbon ()Wevrather et al Ill999l ) and oxygen ions ([Stoll et al 1 11995) 
have been gathered over wide energy (and LET) ranges of these particles, corresponding 
to penetration depths of approx. 45 /xm to 140 mm for carbon ions and 30 /im to 200 
mm for oxygen ions, covering the majority of clinically relevant ranges. Figures 0] and El 
illustrate that the simplified probabilistic model enables to represent cell survival over 
such wide ranges in a systematic manner by using a low number of parameters only (the 
present analysis involves 6 parameters, a ,ai,a 2 and b ,bi,b 2 , plus the cross section of 
cell nuclei, a). This is a prerequisite if the model is to be used in treatment planning 
applications. In the present analysis, cell inactivation effects have been studied for 
monoenergetic particles as function of their LET value (corresponding to their energy). 
For treatment planning applications, an adequate physical model, describing the spectra 
of particle energy and LET values as function of penetration depths and including also 
energy-loss straggling, lateral scattering and fragmentation processes, has to be used 
together with the model of biological effects presented here. Such a phys ical model 
might be based on detailed Monte Carlo calculations (e.g. (IZieglerl 12004)) or semi- 
phenomenologic al parameterization of measured experimental data ([Kramer et al 2000, 



Chu et al 1119931 ). Calculations of survival as function of penetration depth for several 



cell lines, using the radiobiological module presented here and a simplified description 
of physical processes based on SRIM-2003 calculations, have been already performed, 
showing good a greement with measured data; the results will be presented separately 
( Kundratl 12005). To enhance the predictive power of the model, relation between the 
damage induction probabilities a, b and track structure characteristics of different ions 
will be sought, with the aim to help in identifying the optimal ion in different clinical 
situations. Additional systematic experimental data, especially for inactivation induced 
by light ions (helium to nitrogen), would be helpful in solving this issue, too. 



5. Conclusion 



Systematic analysis of published survival data for V79 cells irradiated by light ions has 
been performed using a simplified scheme of the probabilistic two-stage model. The 
probabilities of single ions to induce severe damage to DNA have been derived for 
different ions in dependence on their LET values. In the lower-LET regions, combined 
damages dominate, while survival at higher LET values is governed by single-particle 
damages (Figure [HJ). The present results give quantitative estimates of the increase in 
damage complexity with increasing atomic number Z and LET value. The analysis of 
survival data for CHO-K1 cell line shows then the differences in DNA damage induction 
for different cell lines. 

The derived probabilities of inducing single-track lethal damages have been 
compared with the Monte Carlo estimates of the yields of complex DNA lesions. Good 
agreement between these two methodologically different approaches was demonstrated, 
indicating the biophysical interpretation of lethal events and the importance of damage 
complexity with respect to its biological consequences. 
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The results of this work will be used in proposing detailed biology-oriented 
approaches for hadrontherapy treatment planning. The results may also contribute to 
understanding the differences between the mechanisms of biological effects of different 
ions. 
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